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ABSTRACT 

The random-bond XY spin glass with ferromagnetic next-nearest-neighbour interactions 
is studied on a square lattice by Monte Carlo simulations. We find strong evidence for 
a finite-temperature spin glass transition at Tc ^ 1.1. We also give estimates for the 
spin glass critical exponents for different values of the strength of the nearest-neighbour 
interaction. Our results are consistent with universal behaviour. 



2 



Over the past two decades there has been considerable interest in spin glasses. There 
now exists substantial numerical evidence that the spin glass transition occurs at zero 
temperature in two dimensions (2d) for both Ising [1-5] and vector [6-9] spin glasses with 
nearest-neighbour interactions. 

Spin glass order is widely believed to occur at a finite temperature for the Ising spin glass 
in three dimensions [3,4]. Hence, the lower critical dimension of the Ising spin glass is 
expected to be between 2 and 3 [10,11]. 

The two-dimensional Ising spin glass has been re-investigated in recent years [12,13]. Shi- 
rakura and Matsubara [12] have presented numerical evidence to suggest that the asym- 
metric random bond {+J and — 0.8J) Ising spin glass in 2d has a finite-temperature spin 
glass phase transition. 

Further evidence for a non-zero transition for the Ising spin glass in 2d came re- 
cently from a Monte Carlo study by Lemke and Campbell [13] who modified the Edwards- 
Anderson model to include ferromagnetic interactions between next- nearest-neighbours. 
Lemke and Campbell [13] found that the modified 2d Edwards- Anderson Ising spin glass 
exhibits behaviour very similar to that seen in conventional Ising spin glasses in 3d [4] . 
For vector spin glasses, such as XY spin glasses, the situation is somewhat more contro- 
versial and complicated because of the presence of chirality. The numerical evidence [6-9] 
clearly points to a zero-temperature transition in the nearest-neighbour XY spin glass in 
2d with both Gaussian and random (± J) bond distributions. The chiral-glass transition is 
also believed to occur at zero-temperature [8,14,15]. However, it has been argued that the 
chiral and phase variables decouple on long length scales. As a consequence, the values of 
the chiral- and spin-glass correlation- length exponents are significantly different [8,15]. 

Although most work would seem to point towards a zero-temperature phase transition 
also in three dimensions [6,7,16], some existing Monte Carlo data can be fitted equally well 
assuming a finite temperature transition [7] . 

Very recently, Maucourt and Grempel [17] have used a domain wall renormalisation 
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group method to suggest that the lower critical dimension of the XY spin glass model 
is very close to three, even possibly three itself. The results, however, point to a finite- 
temperature chiral-glass transition in 3d [17]. 

Given the recent developments for Ising spin glasses, in this Rapid Communication we 
re-visit the XY spin glass in 2d. Here we concentrate on the spin glass transition. We 
shall present numerical evidence that including next-nearest-neighbour ferromagnetic in- 
teractions in a random-bond XY spin glass in 2d induces a finite-temperature spin glass 
transition. 

The Hamiltonian for the modified model [13] we study is given by 

H = — JjjSj.Sj — y ^ Sk-Sm, (1) 

<i,3> «k,m>> 

where Sj are planar spins on a square lattice of size {L < 12) with periodic boundary 
conditions. Whereas the first summation < i,j > runs over all nearest-neighbour pairs 
only, the second summation << k,m » denotes sums over all second nearest- neighbour 
pairs a distance y/2 apart. The interactions Jy are quenched independent random variables 
chosen from the following binary distribution: 

PiJij) = l[SiJij-X) + SiJ,j + X)]. (2) 

In our simulations we consider < A < 1. Clearly, for A = the lattice decouples into two 
independent inter-penetrating square sub-lattices. Each sub-lattice is a pure XY model 
which will undergo a Kosterlitz-Thouless [18] phase transition at T ^ 0.89 [19]. 
In our Monte Carlo simulations we study the dimensionless Binder parameter gsG given 

by [4] 

(4) 

»- = 3-2jfc. (3) 

where q^sG^ ^^e spin glass order parameter, and are defined by 

1S^=^J2^<S^■^J>T]J (4) 



and 

4g=]^ E [<S..S,-S„.Sp>|]j (5) 

i,j,n,p 

and here < ... >t denotes a thermal average, indicates an average over disorder and 
N = LxL. 

(2) 

The spin glass susceptibility, xsGi is related to q^Q by 

XSG = Nqfl. (6) 

We analyse our results according to standard finite size scaling [4]. Near a transition 
temperature Tc the Binder parameter is expected to scale as 

9SG{L,T)=gsG{L'/''{T-T,)), (7) 

where g^Q and v are the scaling function and the correlation length exponent, respectively. 
The scaling form of the Binder parameter can be used to determine the value of Tg as 
gsG{L,Tc) = ^5g(0) is independent of the system size L. 
The scaling form of xsG is given by 

XsG{L,T) = L'-^XsGiL'/''{T-T,)), (8) 

where XsG i^ scaling function and rj is the exponent describing the power-law decay of 
correlations at T^. The value of ry can be determined from 

XSGiL,Te)^L^-'^XsG{0)- (9) 

We now turn to our computer simulations and discuss the results. We use Metropolis 
dynamics and sequential updating in our Monte Carlo simulations. The method of Bhatt 
and Young [4] is used to ensure that thermal equilibrium is achieved in the simulations. 

For each value of L (4 < L < 12) and A (0.0 < A < 1.0) we averaged over 500 pairs 
of samples. As most of the simulations were performed over a relatively high temperature 
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range (0.7 < T < 1.4), we did not encounter any serious equilibration problems. In fact, we 
were able to achieve equilibrium within 6400 Monte Carlo steps for the largest lattice at the 
lowest temperature considered. Nevertheless, the simulations presented in this work took 
approximately 3 months of CPU time distributed over 10 Silicon Graphics workstations. 

Plots of the Binder parameter for various values of A are shown in figures l(a)-(e). 
The statistical error-bars are in most cases smaller than the size of the data points. Each 
figure shows the data against the temperature for the four different values of L (L = 4, 6, 8 
and 12). The curves for A = 1.0 (figure 1(a)), A = 0.7 (figure 1(b)), A = 0.5 (figure 1(c)) 
and A = 0.3 (figure 1(d)) clearly intersect and splay out below the intersection point. This 
strongly suggests a spin glass transition. Both the transition temperature, Tc(A), and the 
value of the Binder parameter, gsciTc), would appear to depend only marginly, if at all, 
on A. Although there is some uncertainty in the intersection points for all of the values 
of A, our data are consistent with a universal (i.e. A- independent) transition temperature. 
This is to be contrasted with the non-universal behaviour found by Lemke and Campbell 
[13] in the Ising version of the model. The values of Tc(A) are summarised in Table 1 which 
also contains our estimates for the critical exponents (see below). 

The curves in figure 1(e), which shows the plot for A = 0.0, coalesce at around T ^ 0.9 
and then remain together for all lower temperatures. (For this particular value of A we 
went down to a lower temperature, T = 0.5.) This is, of course, what one expects to occur 
at the Kosterlitz-Thouless [18] transition. 

We next extract the critical exponents from the data. For illustrative purposes, we show 
the scaling plots for A = 0.5 only. Assuming that Tc(A = 0.5) = 1.1, in figure 2 we give 
a log-log plot of XSG against L. The slope of the straight line yields a rather large value 
of ?7 = 0.96 ± 0.05. The values we obtained for r] for the other cases of A can be found in 
Table 1. Although there appears to be a lot of scatter in the values for r]{X), given the 
uncertainties in Tc(A), our estimates are not incompatible with a universal value. 

By setting the critical temperature to Tc(A = 0.5) = 1.1 and considering u as an 
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adjustable parameter, we estimate the correlation-length exponent to be v{\ = 0.5) = 
0.75 ± 0.1. The error quoted here is just an estimate of the range of values for which 
the data scale well. In figure 3 we show a scaling plot of the Binder parameter against 
(T - 1.1)lV0-75_ We note ^j^^t this value of v is very similar to that found recently [20] 
for the four- dimensional XY spin glass. Our estimates for the correlation-length exponent 
for the other values of A are given in Table 1. Once again, our results point to a universal 
value of v. 

Finally, as a consistency check, in figure 4 we display for A = 0.5 a scaling plot of 
Xsg/L'^~'^ against (T-Tc)L^/^ with Tc = 1.1, r] = 0.96 and u = 0.75. The large error-bars 
on this plot are a consequence of both the statistical errors in XSG ^ind the uncertainty in 
the value of ij. Nevertheless, all of the data appear to scale reasonably well. 
To conclude, we have presented data from Monte Carlo simulations of a short-range XY 
spin glass in 2d with ferromagnetic next- nearest-neighbour interactions. We find strong 
evidence for a finite-temperature transition. We have used finite-size scaling to estimate 
the spin glass transition temperature and the critical exponents as functions of the strength 
of the nearest-neighbour interaction. Our results are consistent with universal behaviour. 
Furthermore, the value of the critical temperature and the correlation-length exponent are 
very similar to those found recently for the four-dimensional XY spin glass. 
We would like to thank Matthew Birkin for looking after the Silicon Graphics workstations. 
KJH would like to acknowledge the University of Derby for a Research Studentship. 
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Table 1 

Estimates of the critical temperature, Tg, and the critical exponents, rj and ly, for various 
different values of A. 



A 


Tc 


7] 


V 


0.3 


1.10 ±0.05 


0.85 ± 0.06 


0.85 ± 0.05 


0.5 


1.10 ±0.05 


0.96 ±0.05 


0.75 ±0.1 


0.7 


1.05 ±0.10 


1.08 ±0.06 


0.75 ± 0.05 


1.0 


0.90 ±0.15 


0.83 ± 0.05 


0.90 ±0.1 
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F 

IGURE CAPTIONS 
Figure 1(a) 

A plot of the Binder parameter, qsg^ against the temperature for L = 4, 6, 8 and 12 with 
A = 1.0. The hues are just guides to the eye. 

Figure 1(b) 

The Binder parameter against the temperature for A = 0.7. The lines are just guides to 
the eye. 

Figure 1(c) 

The Binder parameter against the temperature for A = 0.5. The lines are just guides to 
the eye. 

Figure 1(d) 

The Binder parameter against the temperature for A = 0.3. The lines are just guides to 
the eye. 

Figure 1(e) 

A plot of gsG against the temperature. The lines, which are just to guide the eye, clearly 
meet and then remain together for lower temperatures. 

Figure 2 

A log-log plot of the spin glass susceptibility against the linear dimension of the lattice for 
A = 0.5. The slope of the straight line gives us a value of 77 = 0.96 ± 0.05. 



Figure 3 

A scaling plot of gsG against (T 



- Tc)Li/^ with Tc = 1.1 and v = 0.75 (here A = 0.5). 
9 



Figure 4 

A scaling plot of xsg/L'^~'' versus (T-Tc)Li/^ for A = 0.5 assuming that = 1.1, u = 0.75 
and r] = 0.96. The solid line is just to guide the eye. 
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Figure 1 (b) 




Figure 1 (c) 




Figure 1 (d) 




Figure 1 (e) 
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